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I. INTRODUCTION 



Recently charge symmetry violation (CSV) in the parton distributions of the nucleon 
has attracted great interest. It has been generally believed that charge symmetry (CS) was 
highly respected in the nucleon system. Most low energy experiments have shown that CS 
is satisfied within about 1% in reaction amplitudes and most high energy tests are also 
consistent with the charge symmetry, although generally with less precision than the low 
energy tests (for a recent review see and references therein). However some theoretical 
calculations have suggested that the CSV in the valence quark distributions may be as large 
as 5% — 10% |0-§[ which is rather large compared with the low-energy results. There have 
been proposed a number of experiments in which CSV may be observed 

A serious challenge to CS has come from the comparison of the F2 structure functions 



measured in charged and uncharged lepton deep inelastic scattering |TT] , p!0| performed by 
Boros, Londergan and Thomas |T^. A significantly larger CSV than the expectations of 
both theory and other experiments |jl3| was found for the s and s distributions in the 
low-x region (x < 0.1). Any unexpected large CSV will greatly affect our understanding 
of non-perturbative dynamics and hadronic structure 0, and also the extraction of sin^Oy/ 
from neutrino scattering More recently, Boros, Steffens, Londergan and Thomas |jl4 



performed a similar analysis to Ref. [0 with improved corrections for nuclear shadowing 
and the charm threshold in the neutrino data. They found that the data (including the low- 
X region where a large discrepancy was found in [Q) are consistent with charge symmetry 
within experimental errors and the present uncertainty in the strange quark distribution of 
the nucleon. 

Thus suggestions of any large CSV in the parton distributions of the nucleon result from 
theoretical calculations Most of these theoretical calculations are performed using a 

quark model, such as the MIT bag model @], Los Alamos Potential Model or a model 
independent version f^J^. The quark model calculations are based on a quark-spectator 
(quark-diquark) picture of the nucleon in deep inelastic reactions, which is questionable in 
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the low-x region. Hence the quark model predictions for the low-x behaviour of CSV are not 
very rehable. For the sea quark content, quark model calculations involve spectator states 
containing four quarks, which is also unreliable in the low-x region as the mass parameter 
for these four-quark (or three-quark one-antiquark) states is less well-determined than that 
for the diquark states. CSV has also been estimated using the light-cone baryon-meson 
fluctuation model [Q. However, this calculation is mostly qualitative - the quantitative 
calculations are highly dependant on the model parameters. More theoretical study on CSV 
from a different point of view than the quark model will be worthwhile. 

In this paper, we point out that the CSV in both the valence quark and sea quark has 
the same non-perturbative source, and that the meson cloud model (MCM) - which has 
been successfully applied to the study of the sea quark content of nucleon (including SU(2) 
flavour asymmetry |]T^|-[T^, and s and s content of the nucleon [rB,2U|,^) - can provide a 
natural explanation of CSV in the valence quark and sea quark distributions of the nucleon. 
We shall make an alternative investigation of CSV in the parton distributions of the nucleon 
by using the meson cloud model instead of the more commonly used quark model. Our 
calculations for the CSV are significantly different from the quark model predictions. 



II. CHARGE SYMMETRY VIOLATION IN THE MESON CLOUD MODEL 

Charge symmetry results from the 180°-rotation invariance of the strong Hamiltonian 
about the 2-axis in isospace 0. At the quark level charge symmetry implies the invariance 
of a system under the interchange of up and down quarks. For the valence and sea parton 
distributions, this results in the following relations: 

ul{x) = d:{x), rf^(x)=<(x), (1) 
uP{x) = r{x), dP{x) = m"(x), (2) 
sP{x) = s"(a;), s^ix) = s"(x). (3) 

The charge symmetry violation in the parton distributions of the nucleon can be 'measured' 
via the quantities: 
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5d,ix) = - <(x), 6u,ix) = - d^ix), (4) 

Sd{x) = dP{x)-u''{x), 6u{x) =u%x) -d'^ix), (5) 
6s{x) = sP{x) -s'^ix), 6s{x) = sP{x)-s''{x). (6) 

Before entering detailed calculation, it is helpful to break down the parton distribution 
in the nucleon into three parts: the parton distribution in the bare nucleon, the perturbative 
contribution, and the non-perturbative contribution i.e. 

(jp — (jv A_ I ^p o.p — n.p I ..p I ..p C7^ 

" "bare "per "non' " "bare "per "non' \' ) 

" "per ^ "non' " "per ^ "non' V~') 
^ ■'per ^ ''non' ''per ^ ''non- ^^z* 

Similar relations exist for the parton distribution of the neutron. We expect that the bare 
part obeys the charge symmetry 

^bare ~ "^bare' "^bare ~ '^bare" (10) 

The perturbative sea is produced in a very short time via gluon splitting, thus we expect 
the perturbative sea to also be SU(2) flavour symmetric. 
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From the bare parton distribution being charge symmetric [Eq. (|TOD] and the perturbative 
sea being quark-antiquark symmetric [Eq. (p!^] and charge symmetric [Eqs. (p!3| ) and (0)] 
we have the CS violating valence distributions 



6dy — {d'^on ^non) i'^non '^non) ^ (IS) 
SUy = ("U^Qjj — Mnon) ^ ('^non ~ '^non)- i^'^) 

Using the charge symmetry of the perturbative sea [Eqs. (|13|)-(|1^)] we can obtain the CS 
violating sea distributions 

Sd = - Su = ul^^ - (18) 

Ss = - s" Ss = s'P - (19) 

Thus the charge symmetry violation in both the valence and the sea distributions has a 
non-perturbative origin. 

The meson cloud model (MCM) is a model of the non-perturbative contribution to the 
quark distributions of the nucleon. It can provide natural explanations of the flavour asym- 
metry in the nucleon sea [jl5|-0 and quark-antiquark asymmetry [l^,|20|j2l| in the nucleon. 



The essential point of the MCM is that the nucleon can fluctuate into different baryon-meson 
Fock states, 

|iV)physicai = ^|iV)bare + E E / d''i^± kl) \B\y, k^); M^'(l - y, -k^)) (20) 

where Z is the wave function renormalization constant, 'Psifiy, kj_) is the wave function of 
the Fock state containing a baryon {B) with longitudinal momentum fraction y, transverse 
momentum k_L, and helicity A, and a meson (M) with momentum fraction 1 — y, transverse 
momentum — k^, and helicity A'. The model assumes that the lifetime of a virtual baryon- 
meson Fock state is much larger than the interaction time in the deep inelastic or Drell-Yan 
process, thus the contribution from the virtual baryon-meson Fock states to the quark and 
anti-quark distributions of the nucleon can be written as convolutions 



qnon{x) = E 



BM 

Q'non('^) 



fBM/N{y)q^{-)+ / -fMB/N{l-y)q'\-) 
X y y Jx y y 



(21) 



E t-fMBMl-y)q'\-), (22) 

BM ■^'^ y y 

where fBM/N{y) = fMB/Ni)- — y) is fluctuation function which gives the probability for the 
nucleon fluctuating into a virtual BM state 



oo 
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fBM/N{y) = J^ dk_^ (l)BM{y:k_J . (23) 

As the proton and neutron form an SU(2) isospin doublet, the baryons and mesons in their 
respective virtual Fock states differ only in their carried charge. If we neglect the mass 
differences among these baryons and mesons the fluctuation functions for the proton and 
neutron will be the same. Thus the contributions to the parton distribution of the nucleon 
from these fluctuations will be isospin symmetric (charge symmetric). As is well known, the 
electromagnetic interaction induces mass differences among these baryons and mesons. If we 
take into account these mass differences, the probabilities for the corresponding fluctuations 
of proton and neutron will be different and thus the contributions to the parton distributions 
of the proton and neutron will be different, which results in CSV in the parton distributions 
of the nucleon. Thus the MCM can provide a natural explanation of CSV in the parton 
distributions of the nucleon. Although it is argued from the quark model calculations that 
the electromagnetic effect does not play a signiflcant role in the calculation of CSV in the 
parton distributions, it is worthwhile to study this effect using a different theoretical picture. 

For the CSV in the up and down quark distributions, we consider the fluctuations N — > 
Nn and N — > Att, but neglect the other fluctuations such as N ^ N (A)p, N ^ Nrj{u;) 
and Ar] {u) due to the higher masses of the involved mesons. Thus the fluctuations we 
consider include: 

p{uud) n{udd) + tt "*"(«(!), n{udd) p{uud) + T:~{ud), 
p{uud) — >• A^{udd) + 7i~^{ud), n{udd) — >• A'^{uud) + 7i~{ud), 
p{uud) p{uud) + Tx^{[uu + dd]/\/2), n{udd) n{udd) + 'k^{[uu + dd]/\/2), 
p{uud) A'^{uud) + 7r°([ii'u + dd]/V2), n{udd) A^{udd) + 'k^{[uu + dd]/^/2), 

p{uud) A++(tititi) + 'K~{ud), n{udd) A~{ddd) + 7r+(ti(i). (24) 

For the CSV in the strange and anti-strange quark distributions we consider the following 
fluctuations 

p{uud) A{uds) + K~^{us), n{udd) A{uds) + K°{ds), 



p{uud) S°(mc/s) + K+{us), n{udd) TP{uds) + K^{ds), 
piuud) — > Ti^iuus) + K^{ds), n{udd) T,^{dds) + K^{us). 



(25) 



Since we consider the CSV between the proton and neutron, we should not neglect the mass 
difference between the proton and neutron, nip — rrin = — l.SMeV, and those among the 
baryon and meson multiplets ||22|| , 



mA- — m^o = m^o — m/^+ = — m^^ 

m-s- - mso = 4.8 MeV, 

m-^o — ms+ = 3.3 MeV, 

m^ri — rrij^o = 4.6 MeV, 
rriKO — mx+ = 4.0 MeV. 



1.3 MeV, 



(26) 



The probabilities of various fluctuations can be calculated using the effective Lagrangian 
and time-ordered perturbation theory in the infinite momentum frame. For the fluctuations 



— > A^TT, KK and — * UK, the fluctuation functions can be expressed as |[16| , p!8 



fBM/N{y) = Ci'- 



167r2 Jo y{l - y) {m% - m%j^jy 



y 



(27) 



where y is the longitudinal momentum fraction of the baryon B, Qnbm is the effective 
coupling constant, ttl'^bm is the invariant mass squared of the BM Fock state. 



m 



BM 



y 1 - 2/ 

and Gbm is the phenomenological form factor, for which we adopt the exponential form 

'm% -ml^j{y, kl) 



{2i 



GsMiy.k^) = exp 



2A2 



(29) 



A is a cut-off parameter which can be taken as A = l.OSGeV for all fluctuations involving 



octet baryons and pseudoscalar or vector mesons [T^. The effective coupling constants are 
taken to be Qnntt = 13.07 |[T6|jr8[| , Qnak = 13.12 and gNT:K = 6.82 [^^ . The coefficient 
Ci in Eq. (p7|) comes from the Clebsch-Gordan coefficients for the fluctuations of different 
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isospin multiplets. Ci = I for /p^o/p, /„^o/„, /A/r+/p, /axo/™, h^K+ip and h-K^/n^ and 
Ci = 2 for fnn+/p, fpn-/n, fs+RO/p, and f-£-K+/n- For the fluctuation ^ Att we have 



^'^''^ ^ 167r2 Jo y{l — y) (m^ — m\^Y Qm\y^ 



where g = 11.8 GeV ^ |jT^,[T8|, C2 = I for f/\++-n-/p and fA-TT+/n, C2 = 2/3 for /a+ttO/p and 
/AOvrO/n, and C2 = 1/3 for fAOn+/p and fA+n-/p- We adopt the exponential form [Eq. (p9|)] 



for the form factor and the cut-off parameter was taken to be A = 0.98 GeV |]I8| . 

In the meson cloud model the non-perturbative contribution to the quark and the anti- 
quark distributions in the nucleon sea come from the quarks and anti-quarks of the baryons 
(A^, A, S) and mesons (vr, K) in the virtual baryon-meson Fock states. So we need the parton 
distributions in the involved baryons and mesons as input. For the parton distribution in 
the pion, we employ the parameterization given by Gliick, Reya, and Stratmann (GRS98) 



24|] and we neglect the sea content in the meson, that is, 



(J'^+ = = u""' = (F~ = ^v"", (31) 

= = dr' = (T' = (32) 

v^ix, /x^Lo) = 1.052x-°-^9^(l + 0.357v^) (1 - a;)°■3'^^ (33) 

at scale /i^Lo ~ 0-34 GeV^. For the s distribution in the and K'^ we use the GRS98 
parameterization 



s''\x,^'^^o) = ^"""(a^^/^NLo) = [1 - 0.540(1 - xy>-''\ t;-(x,/x^Lo) (34) 

at scale /i^Lo ~ 0-34 GeV^. For the quark distributions in the bare baryons, we first use the 
up and down quark distributions in the proton given by Gliick, Reya, and Vogt (GRV98) 



d'{x, /x^Lo) = 0.624(1 - x)u^{x, /i^Lo) (35) 
u^{x, /^NLo) = 0.632 x-^-^\l - x)2-°^(l + 18.2x), (36) 
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at scale /i^Lo ~ 0-^0 GeV^, then relate these to the distributions in the other baryons via 
the relations 



u 



A+^ 



U 



AO 



uP + dP, =uP + dP, 

T- 1 



(37) 
(38) 
(39) 



We evolve the distributions to the scale = A GeV^ using the program of Miyama and Ku- 
mano 



26| in which the evolution equation is solved numerically using a brute-force method. 
We found that at = 4 GeV^ the parton distributions we need (f '^(x, Q^), s^°{x,Q'^) 



uP{x, Q"^) and (i^(x, Q"^)) can be parametrized using the following form 

q{x,Q'^) = ax'' (1 - xY {l + d^/x + ex) 



(40) 



with the parameters given in Table 1. We estimate the uncertainty in solving the evolution 
equations numerically and parametrizating the parton distribution in the form of Eq. ( ^OD 
to be about 2% in the x-region which we are interested in ie x > 10"^. 

The final expressions for the CSV in the valence parton distributions are given by: 



x6d^ 



^ X I V 1 X 

dy- { fnn+/p{y) - fpn-/n{y) M^(-) 
y ' ' ^ y 

1 1 X 

/n7r+/p(l - y) - /p7r-/n(l " y) T^^'"!") 
1 X 

fpnO/piy) - fnnO/n{y)\ <P{-) 

1 X 

fAOn+/p{y) - fA+n-/n{y) M^(-) 
' i y 

/A07r+/p(l -y) - /A+7r-/n(l - 1/) 



+ 
+ 

+ 
+ 



2 y 



fA+wo/piy) - fAo-no/niy) d^i 



,x. 



/A++7r-/p(l -y) - /A-7r+/n(l - Z/) 



2 y 



(41) 



^ X \ r ~\ X 

dy- { fn^+/p{y) - fp7T-/n{y) d^i-) 
y ^ y 



1 „ ,x. 



n7r+/p\ 



1 -y) - /p7r-/n(l -y) o^^(- 



+ 
+ 

+ 
+ 
+ 



fpnO/p{y) - fnnO/n{y) 



,x , 



fAo^+/piy) - fA+n-fniy) dP{-) 

' i y 

/A07r+/p(l -y) - /A+7r-/n(l - ?/) 



2"''? 



fA+nO/piy) - fAOTrO/niy) 



,X , 



fA++n-/p{y) - fA-7T+/n{y)] + dP) 

-1 1 X ] 

/A++7r-/p(l -y)- /A-7r+/n(l " y)\ if'' \ ' 



(42) 



For CSV in the sea we obtain 



x5d 



^ X ( 1 r 

^ dy- |/„^+/p(l -y) - /p7r-/n(l ~y)+ 2 -y) - fnnO/ni^ " 1/) 



+ /A07r+/p(l - y) - /A+7r-/n(l - V) 
1 



+ 2 



1 ^,X, 



(43) 



X r 1 



X^M = dy- |- [/p^O/p(l -y)- fnnO/n{i " Z/) 
+ /A++7r-/p(l -y) - /A-7r+/n(l " 



+ 2 



/A+7r0/p(l -y) - /A07r0/n(l " ?/) 



(44) 



x5s = dy- [fAK+/p{y) - fAKO/n{y) + hoK+/p{y) - hoKO/n{y) 



+fT.+KO/p{y) - f^-K+/n 



1 T 

2 1/ 



(45) 



f X ( 

XSS = dy- |/Ai^+/p(l -y) - fARO/ni^ -y) + hoK+/p{^ -y)- hoKO/ni^ " Z/) 

+h+KVp{i -y)~ h-K+fnii - y)} Is^'i^) (46) 



III. RESULT AND DISCUSSION 



From Eqs. (^if) - (l46| ) we can see explicitly that the differences among various ffuctuation 
functions such as fnn+/p and fpir-fn result in the CSV in the parton distributions of the 
nucleon. We plot these differences in Figs. 1-3. It can be seen that the difference in the 



fluctuation — > iVvr {fnn+/p — fpn-/n, •••) is much larger than that in the fluctuation 
— > Att (/A++7r-/p ~ /A+7r+ /n, ••■) and the latter is much bigger than that in the fluctuation 
— > AK (/att+Zp — fARO/n, •••)■ Thus the CSV in the valence and sea up and down quarks 
should be the same order, and both larger than the CSV in the s and s distributions. The 
difference fmr+fp — fp-K-fn is much larger than fpT,o/p — fmrO/n, thus the CSV in the sea of 
the minority quark flavor {6d) will be much larger than that of the majority quark flavor 
{6u) due to the absence of {fmr+fp — fp-K~/n) term in 5u. The probabilities of the various 
fluctuations can be obtained by integrating the corresponding fluctuation functions. We flnd 
the probabilities of the dominant fluctuations to be 

P{j) n7r+) = 0.202 P{n pn') = 0.205, (47) 
P{p A++7r-) = 0.0481 P{n A-n+) = 0.0475, (48) 
P{p AK+) = 0.0127 P{n AK°) = 0.0125, (49) 

that is there is about a 1% excess of fluctuations n pn^ over p mi^ and p A++7r~ 
over n A^7r+, and about 2% excess of p KK^ over n —>■ AK^. 

We present our results for the CSV in the valence quark sector {x6d^ and xSu^) in 
Fig. 4. We flnd that x6d^ and x6uy have similar shape and both are negative, which is 
quite different from the quark model prediction of xSd^ being positive for most values of 
X Furthermore, our numerical results are about 10% of the quark model estimation 

[0-01 . It has been argued that although the absolute values of 6d^ and 6uy are small, the 
ratio -Rmin = Sd^/d^ may be much larger than the ratio -Rmaj = ^Uy/u^ in the large-x region 
since the dP{x) / uP{x) ^ 1/2 as x — 0, and values as large as 5% ~ 10% have been 
obtained for the ratio Sd^/d^. No such large-x enhancement appears in our calculation for 
both ratios. We flnd that the ratio 5dy/d^ exhibits a maximum about 0.2% at x = 0.1 while 
the ratio du^/u^ diverges as x but is smaller than 0.3% in the region of x > 0.02. The 
numerical results for the CSV in the sea quark {x5d, x6u, x6s and x6s) are given Fig. 5. 
We flnd that x6d has the largest CSV and that xSd^ and xSu^ are of similar magnitude, 
which is consistent with our expectation from the analysis of the fluctuation functions. Our 
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prediction for the x^J being negative is opposite to the positive theoretical prediction in 0. 
Our calculation for the low-x behaviours of x5d, x6u, x6s and x6s are also quite different 
from the quark model prediction - the quark model predicts that these quantities diverge 
as X ^ while our calculations show all these CSV distributions go to as a; 0. We 
did not find any significant large CSV in the sea quark distribution of the nucleon, which is 
consistent with the most recent phenomenological analysis |M . We would like to emphasise 



that instead of the quark model we adopt a totally different model, the meson cloud model, 
to calculate the charge symmetry violation in the parton distributions of the nucleon. The 
quark model calculation in the small-x region is not very reliable since the quark-diquark 
picture that is employed breaks down in this region. The meson cloud model is suitable 
in the study of the CSV in the parton distribution of the nucleon since it has the same 
non-perturbative origin as the d/u asymmetry in the proton. 

IV. CRITICISM OF QUARK MODEL CALCULATIONS 

We have already mentioned a few of the difficulties with CSV calculations using quark 
models. A recent paper by Benesh and Londergan |^ attempted to avoid any quark model 
specifics and relate possible CSV in the valence quark distributions to the measured valence 
distributions. Starting from the parton model expression for a quark distribution 

q{x) = p+ ^ 6{p+{l -x)- p+) I {n\^{0) \p) f (50) 

n 

where the intermediate state \n) has 4-momentum pn and the + components of momenta are 
defined by = + k^, and then making the assumption that the intermediate state can 
be modelled by a diquark system with definite mass Md, Benesh and Londergan investigate 



the consequences for CSV of varying M^. Following the Adelaide group p8[, we can attempt 
to determine the dependence of the quark distribution on Md. Assuming that the modulus 
squared of the wavefunction for the struck quark in the nucleon is symmetric about the 
z-axis, we can use the delta function to perform the integration over transverse diquark 
momenta 
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/roo 
(ip„5(p+(l - x) = 2tt dpnPn (51) 

Pmin 



where 



Pn 



M\l - xf - Ml 



2M(1 - x) 



(52) 



PT = 2M(1 - x)^Mj + p2 - M\l - x)' - Ml (53) 
and M is the nucleoli mass. Therefore we obtain the quark distribution in the form 



= / dpngiPn) (54) 

where g{k) only depends on the magnitude of the 3-momentum (this is not true in the case of 
spin dependent quark distributions), and we have reminded ourselves that Pmin is a function 
of X and Mj_. Thus all the M^ dependence of the quark distribution is in the lower limit of 
the integral. Use of the fundamental theorem of calculus then gives 

-/^r- (55) 



dMd dx dMd ' dx 

2Md{l-x) dq{x) 



(56) 



M2(l-x)2 + Mj dx 

This expression is similar to that of Benesh and Londergan |^ , except that in their case the 
d/dx operator acts on the product of the kinematic factor 2Md{l — x)/{M'^{l — x)^ + Mj) 
and the original quark distribution. The derivation of reference differs from ours in that 
they make a variation in M^ under the integral in equation (^0]), then evaluate the integral 
over Pn by ignoring any transverse momenta of the diquark. However in our expression, all 
transverse momenta have been properly integrated over (in the parton model the transverse 
momentum of the struck quark vanishes). 



Benesh and Londergan then use the idea of Close and Thomas |2^ that the quark model 
SU{4) spin-isospin symmetry is broken by the color hyperfine interaction. The hyperfine 
interaction leads to a splitting in the masses of the spin-0 and spin-1 diquark states and 
hence to a difference between the up and down valence distributions: 

3 1 

d,{x) = qlix) (57) 
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where the superscripts s, t refer to singlet and triplet diquark states respectively. If the 
N — A mass splitting is also caused by the color hyperfine interaction, then the shifts in the 
singlet and triplet diquark masses are found to be 



1 

r 



+50MeV. 



(58) 



By now expanding M^) and M^) in Taylor series to first order in 5Md about the 
symmetry point = M°, q^{x,M^) = g*(x,M°) Benesh and Londergan obtain the shift 
in the triplet quark distribution 



Sh/qlix) = i(2(i^(x) - Uy{x)). 
o 



(59) 



Now as the CSV at the quark level comes from quark mass and electromagnetic effects, both 
of which are iso- vector, the only mass shift is in M^, and to first order the shift in the triplet 
quark distribution will be proportional to that from the hyperfine interaction 

5csvM\ 2d^{x) - Uy{x) 



Scsvqvix) 



(60) 



The main difficulty with this argument is that it is entirely based on first order shifts 
in the quark distributions. However the second order terms can be estimated, and they are 
of similar magnitude to the first order terms. Expanding the quark distributions to second 
order in SM^, about the symmetry point we have 



q,{x, MO + 6Md) = g.(x, M°) + 6M, 



dqy(x,Md) 
dMd 



2 d'^qy{x,Md 
d{Mdf 



where the partial derivatives on the right hand side can be evaluated using 
gives for the hyperfine shift in the triplet quark distribution 



. This then 



1 



5hfql{x) = --5ufql{x)+2{5M^ 



t\2 



d\y{x,Md) 



diMd] 



1/ , / N / NN 3 , ^ , ,9 (7„(a;, Mrf) 
-(2dy{x)-Uy{x)) + -{6M\f ^''^ 



diMd) 



Ma=MO 



(62) 



As an estimate of the second order term we can approximate qy{x,M^) by d^x or {uy{x) + 
dy {x))/3 and use one of the well-known parametrizations of the quark valence distributions 
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[^,^. We also take our value for M° to lie in the range (0.65 — 0.85)M, which is the range 
suggested by quark models, though the results are not very sensitive to the value of M° in 
this range. In Fig. 6 we compare the first and second order terms for 6hfql{x). We can see 
that our estimate of the second order term is of comparable magnitude to the first order 
term over most of the z-range. Indeed for low x it is larger, showing the unreliability of 
quark model calculations in this region. For x > 0.2 the second order term is of opposite 
sign to the first order term, which indicates that the first order estimate of the shift in 
the quark distribution is too large in this region. This in turn implies that the estimate 
of CSV induced shifts in the quark distributions in this region are also too large. These 
conclusions are not much influenced by the choice of valence quark parametrization, the value 
chosen for M°, or whether we use Benesh and Londergan's expression for the dependence of 
the quark distribution on rather than equation (|5B|). The reason for these conclusions 
not being greatly influenced by the choice of expression for dq{x) / dM^ is that, with the 
quark distributions used, the highest order derivative term in x always dominates. This is a 
consequence of the divergences in the valence quark distributions near x = 0, qv{x) ~ 
in all cases. 



V. SUMMARY 

Although it has been generally assumed that charge symmetry was highly respected in 
the nucleon system, there have been some phenomenological analysis [P!^,p!^ and theoretical 



calculations [^-|^ about the possible extent of CSV in the parton distributions. Any unex- 
pected large CSV will greatly affect our understanding on the non-perturbative dynamics and 
hadronic structure, and the extraction of sin^Oy/ from neutrino scattering. Up to now most 
theoretical attempts to calculate the CSV in the parton distributions are based on the quark 
model and employ the quark-diquark model. In this paper we point out that CSV in both 
the valence and sea quark distributions of the nucleon can arise from the non-perturbative 
dynamics of the nucleon. We present an alternative analysis of CSV in the parton distribu- 
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tioiis employing the meson cloud model, which has previously been successful in the study 
of the flavour asymmetry and the quark- antiquark asymmetry of the nucleon. In the meson 
cloud model the proton and neutron may fluctuate into hadron-meson Fock states in which 
the hadrons and mesons are in different charged states respectively. As we consider the 
mass differences among these hadrons and mesons, the probabilities of proton and neutron 
fluctuating into the corresponding Fock states will be different. Thus the non-perturbative 
contributions to the valence and sea quarks distributions will be different, which naturally 
leads to the CSV in both the valence and sea distributions of the nucleon. Our predictions 
for the CSV in the valence sector and sea sector are both different from the quark model 
calculations. We also point out the deflciencies of quark model based calculations of CSV 
in the parton distributions. In particular the quark-diquark picture is inadequate at low-x, 
and in the medium-x region the use of a flrst order shift in the parton distributions must 
be questioned, as higher order shifts are of similar magnitude. The coming experimental 
information on the parton distributions of the nucleon and more theoretical studies on this 
issue will examine these calculations. 
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Table 1. Parameters in Eq. g^) at = 4 QeV^. 





d 


u 


c 




6 




1.712 


-0.518 


1.182 


-0.836 


0.972 




0.803 


-0.516 


1.306 


-0.762 


0.957 




1.029 


-0.572 


3.933 


1.550 


6.033 


dP{x, Q2) 


0.615 


-0.575 


5.096 


1.102 


6.773 
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FIGURE CAPTIONS 



Fig. 1. The differences between ffuctuation functions Jntt/p and fNn/n- 

Fig. 2. The differences between ffuctuation functions /att/p and fAn/n- 

Fig. 3. The differences between ffuctuation functions /a(s)x/p and /A(i;)x/n- 

Fig. 4 The charge symmetry violation in the valence quark sector. -Rmin = Sdy/d^ and 

Fig. 5 The charge symmetry violation in the sea quark sector. 

Fig. 6 The first and second order shifts in the triplet quark distribution caused by the 
color hyperfine interaction. The solid curve is the first order shift (2^^,(2;) — Uy{x))/6 
calculated using the parametrizations of reference . The dashed curve is the second 
order shift estimated using q{x) = {d^,{x) + u^,{x))/3, = 0.75M, and a mass shift 
of 50 MeV for the triplet diquark state. 
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